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I. INTRODUCTION 



For many years, calculation of Casimir or quantum vacuum energies were confined to 



systems of very high symmetry, such as parallel plates 



2j. Self-energies of spheres |3| and 



cylinders j4| were also calculated. To calculate the force between curved surfaces, such as 
that in the experimentally realizable spherical lens above a plate, resort had to be made to 
the "proximity force approximation" 5|, which, while exact for zero separation, was subject 
to unknown corrections for finite distances between the surfaces. 

That situation has radically changed. The "multiple scattering" technique, not merely 
increase in computing power, has been responsible for much of the improvement. The 



method, related to the Krein trace formula ral, was used early on to construct the Lifshitz 
formula for the quantum fluctuation force between dielectric plates [7|]. It was the basis of the 
classic work by Balian and Duplantier |S] . But recently it was realized by many people that 
it could be applied to practical calculations of Casimir forces between bodies of essentially 
arbitrary shape and orientation j^, Q, H, 12, 13, 14\- For more complete references and 



background, see Refs. [151 . Il6l| . Complementary developments by Gies and collaborators on 
the worldhne method jl?, isilisl, ^^^1 direct numerical methods 21, 22, 23] should also 



be cited. 

In this paper we summarize the work of the Oklahoma group in applying the multiple 
scattering technique to a variety of problems. The hope is that results will be forthcoming 
that would be of use in the design of nanotechnology. We will merely summarize some 
representative results here; for more complete details and further applications the reader is 
referred to the original previous and forthcoming papers. 



II. MULTIPLE SCATTERING TECHNIQUE 

For simplicity, we first restrict attention to the quantum vacuum forces arising from a 
massless scalar field. The multiple scattering approach starts from the well-known formula 
for the vacuum energy or Casimir energy 

^= — TrlnG^ — TrlnGGo\ (2-1) 
2r 2r ° ' ^ ^ 
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where r is the "infinite" time that the configuration exists, and G (Go) is the Green's 
function, which satisfies the differential equation 



(-92 + V)G = 1, -d'^Go = 1, (2.2) 

in terms of the background potential V. 
Now we define the T-matrix, 

T = S -1 = V{l + GoV)-\ (2.3) 

If the potential has two disjoint parts, \^ = Vi + V2, it is easy to derive 

T = {V, + V2){1 - GoT,){l - GoT^GoT^r'il - G^T^), (2.4) 

T, = V,{l + GoVi)-\ 2 = 1,2. (2.5) 

Thus, we can write the general expression for the interaction between the two bodies (po- 
tentials) : 

E12 = -T^Tr ln(l - G0T1G0T2) = -;^Tr ln(l - ^^^11^2^2), (2.6) 
where = (1 + GoVi)-^Go, z = 1, 2. 



A. Exact Results for Weak Coupling 

In weak coupling, where the potentials are very small, it is possible to derive the exact 



(scalar) interaction between two potentials 2J], either in 2- or 3-dimensions: 



where is the "infinite" length in the translationally-invariant direction, and 

SD: ^.-J-,/(*)(.0M^. (2,7b) 

Consider two plates (ribbons) of finite width L, offset by an amount b, separated by a 
distance a, as shown in Fig. [1], 

^i(rx) = Xi6{y)e{x)d{L-x), (2.8a) 
V^ir'^) = X2S{y' - a)e{x' - b)0{L + b - x'), (2.8b) 



3 



b 



L 



a 



L 



FIG. 1: Two ribbons of finite width, offset from each other. The ribbons have infinite extent out 
of the page. 




FIG. 2: Two ribbons of finite width tilted with respect to each other. 

Equation fl2.7ap gives an explicit result for the energy of interaction between the plates: 
E A1A2 

where 



[-2g{b/a) + g{{L - b)/a) + g{{L + b)/a)] , (2.9) 



g{x) =xtan"^x- ^ ln(l + x^) = -Re{l + ix) \n{l + ix) . (2.10) 



We can consider arbitrary lengths and orientations, in 3 dimensions, for the plates [24 1. 

For example, we can consider tilted plates, as shown in Fig. [2l Explicit interaction 
energies can be given in terms of Ti2, the inverse tangent integral. For fixed distance D 
between the center of masses of the plates, and for Li ^ L, L2 00, d — > —00 (that is, the 
upper plate is completely above the much wider lower plate), and if D > ^, the equilibrium 
position of the upper plate is at = 7r/2. That means there is a torque on the upper plate 
tending to orient it to be perpendicular to the lower plate. 

We can also examine the interaction between rectangular parallel plates of finite area, as 
illustrated in Fig. [31 As the perpendicular separation between the plates a tends to 0, the 
force per area on the plates has the expansion 
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FIG. 3: Two parallel finite rectangular plates. 
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FIG. 4: Coaxial disks 



The leading correction to the Lifshitz formula, Ci, has a simple geometrical interpretation.. 
If the upper plate is completely above the lower plate, the leading correction vanishes, 
Ci = 0. On the other hand, if the plates are of the same size and aligned, the correction is 
geometrical: 

1 Perimeter 

Ci = . (2.12) 



TT Area 

In a similar vein, we can consider coaxial disks, as shown in Fig. HI In this case, the leading 
correction is completely consistent with the result for rectangular plates, that is, the leading 
correction vanishes if the disks are of unequal size, -Ri < -R2, ci = 0, while if the are equal, 
Ri = R21 the correction is again given by Eq. fl2.12p . 

We can summarize the salient features for two thin plates as follows: Two plates of 
finite size experience a lateral force so that they wish to align in the position of maximum 
symmetry. In this symmetric configuration, there is a torque about the center of mass of 
a smaller plate so that it tends to seek perpendicular orientation with respect to the larger 
plate. Finally, the first short- distance correction to the normal Lifshitz force is geometrical. 
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FIG. 5: The Casimir pistol. Is it possible to choose parameters so that the "bullet" is expelled 
from the pistol due to quantum vacuum forces? 



The results are of relevance to the recent discussions of the Casimir pistol 25||, shown in 
Fig.O 



III. SUMMING VAN DER WAALS FORCES 



In the preceding section the focus was on scalar field theory. More relevant is elec- 
trodynamics, and the quantum vacuum forces between conducting and dielectric bodies. 
Weak-coupling in this regime is described by the Casimir-Polder force between molecules 
26 1 . The following discussion is based on Ref. [s^]. (For a summary of very extensive work 



in the nonretarded regime, see Ref. 



28|.) 



The (retarded dispersion) van der Waals (vdW) potential between polarizable molecules 
is given by 

V = — : a = - — —■ (3.1) 



47r r'^ ' " A-nN' 

This allows us to consider in the same vein as in Sec. [Ill the electromagnetic interactions 
between distinct dilute dielectric bodies of arbitrary shape. This vdW potential may be 
directly derived from the action 



W 



-Tr inrr;:^ 



where Vi = Sj — 1 and 



V X V X 1- 



3-|CI|r-r'l 

47r|r — r'l 



-l = (VV-lC')Go(r-r'). 



(3.2) 



(3.3) 
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A. Interaction between e, /i bodies 



More generally, we can consider material bodies characterized by a permittivity e{r) and 
a permeability /^(r), so we have corresponding electric and magnetic potentials 



ye(r) = £(r) - 1, Kn(r) = Mr)-l. 



(3.4) 



Then the trace-log in Eq. (13.21) is ($o 



|V X To) 



Tr InTFo^ = -Tr ln(l - ToVe) - Tr ln(l - ToK. 
-Tr ln(l + *oTe$oTj, 



(3.5) 



in terms of the T-matrix, 



(3.6) 



If we have disjoint electric bodies, the interaction term factorizes just as in Eq. (12. 4p : 
Tr In (1 - ToiVi + V2)) = -Tr ln(l + roTi)-Tr ln(l + roT2)-Tr ln(l-roTiroT2), (3.7) 
so only the latter term contributes to the interaction energy. 



Ei,t = -Trln(l-roTiroT2). 
A similar result holds if one body is electric and the other magnetic. 



(3. 



Ei,t = --Trln(l + $oT^$oT^). 



(3.9) 



Using this, it is easy to show that the Lifshitz energy between parallel dielectric and dia- 
magnetic slabs separated by a distance a is 

1 



dC / d'k 



In (1 - nr'^e-^"^) + In (l - r^rse-^'^'^) 



where 



(3.10) 



(3.11) 



with = fc^ + (^^, k\ = k"^ + e^"^, k[ = ki/s, k\ = k"^ + fi('^, n'^ = K2//i. This means that in 
the perfect reflecting limit, e ^ 00, /i — > 00, we get Boyer's repulsive result I29I, 

7 772 



Boyer 



+ 



8 720a3' 



(3.12) 
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B. Dilute dielectrics 



We now give some exact results for dilute dielectrics, |e — 1| ^ 1. For example, consider 
the force between a semi-infinite slab and a slab of finite cross-sectional area A as shown in 
Fig. El The force between the slabs is 









z 












a 




£1 



FIG. 6: Dielectric slab of cross section A above an infinite dielectric plate. 



23 1 



(£l-l)(£2-l), 



(3.13) 



A 6407r2 

which is the Lifshitz formula for infinite (dilute) slabs. Note that there is no correction due 
to the finite area of the upper slab. 

We can also compute the force between a dilute dielectric sphere and a dilute dielectric 
plate, as illustrated in Fig. [3 The energy of interaction is given by 



£2 



R 



£1 



FIG. 7: Dilute dielectric sphere above a dilute dielectric slab. 



23 , ,,47raV3 1 



6407r2 



(3.14) 
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FIG. 8: Rectangular solid (with dielectric constant e^) of side a, 6, and c (perpendicular to the 
plane, not shown) a distance Z above an infinite plate (with dielectric constant £\) extending to 
z = —oo. The shorter side b makes an angle 9 with respect to the plate. 

which agrees with the proximity force approximation in the short separation limit, R — a = 
(5 < a: 

with an exact correction, intermediate between that for scalar l/2(Dirichlet+Neumann) 
[qI. [io| and electromagnetic perfectly-conducting boundaries 11 1. 

1. Torque between slab and plate 

Figure [8] shows a dilute rectangular solid above an infinite dilute plate, where the shorter 
side makes an angle 6 relative to the plane. Generically, the shorter side wants to align 
with the plate, which is obvious geometrically, since that (for fixed center of mass position) 
minimizes the energy. However, if the slab has square cross section, the equilibrium position 
occurs when a corner is closest to the plate, also obvious geometrically. But if the two sides 
are close enough in length, a nontrivial equilibrium position between these extremes can 
occur. Figure [9] shows the equilibrium angle, for fixed center of mass position. For large 
enough separation, the shorter side wants to face the plate, but for 



a 



^<^°=2V ^^-''^ 

the equilibrium angle increases, until finally at Z = D = V + the slab touches the 
plate at an angle 9 = arctan b/ a, that is, the center of mass is just above the point of contact. 



1.0 1.5 2.0 2.5 3.0 3.5 4.0 

ZjD 



FIG. 9: Equilibrium angle as a function of separation of the body from the plane, for given h/a 
ratios 0.95, 0.9, and 0.7, respectively given by solid, dashed, and dot-dashed lines. 




FIG. 10: Dielectric cylinder within hollow parallel cylindrical cavity, 
about which point there is no torque. 



2. Interaction between parallel cylinders 



Two parallel dilute cylindrical bodies of radius a and h (of large length L) , outside each 
other, with a distance R between their axes, have the interaction energy 



E 

T 



23 a%'' 

(£l-l)(£2-l)- 



1 /a^-fe2 

2 I i? 2 



eovr 



5/2 ■ 



(3.17) 



This result can be analytically continued to the case when one dielectric cylinder is entirely 
inside a hoUowed-out cylinder within an infinite dielectric medium js^ , as shown in Fig. [TOl 
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3. Interaction between spheres 



We close this section by giving the interaction energy between two dilute spheres, of 
radius a and b, respectively, separated by a distance R> a + b: 

23 {e,-l){e,-l)( fl-{^-^y 



E 




19207r R 

^ 4aO 53 r O — 52 -L 



i?2 



(3.18) 



This expression, which is rather ugly, may be verified to yield the proximity force theorem: 

It also, in the limit 6 — > oo, i? — oo with R — b=Z held fixed, reduces to the result for the 
interaction of a sphere with an infinite plate. 

IV. EXACT TEMPERATURE RESULTS 

The scalar Casimir free energy between two weak nonoverlapping potentials Vi(r) and 



V2(r) at temperature T is |31| 



T /" , , s , , ,s , X , ,^ coth 27rT I r — r' I , , 

Et = J idv){dr')V,{r)V2{r') _ ' ■ (4.1) 

From this, we see that the free energy of interaction between a semitransparent plane and 
an arbitrarily curved nonintersecting semitransparent surface S, described by the potentials 

V, = \,6{z), V2 = \2S{z'-z{S)), (4.2) 

has the following form: 

Er = [dS [ dx'-^, (4.3) 

where the area integral is over the curved surface. (The upper limit of the x-integral is 
irrelevant, since it does not contribute to the force between the surfaces.) This is precisely 
what one means by the proximity force approximation (PFA): 

i?PFA = j dSS\\{z{S)), (4.4) 

£\\{z{S)) being the Casimir energy between parallel plates separated by a distance ^(5*). 



This is the theorem proved by Decca et al. 32| for the case of gravitational and Yukawa 
forces. 
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A. Interaction between semitransparent spheres 



The free energy of interaction between two weakly-coupled semitransparent spheres, de- 
scribed by the potentials 

Vi = \i5{r-a), V2 = \25{r' -b), (4.5) 

in terms of local spherical coordinates, the centers of which are separated by a distance 
R > a + b, is 

+ f{27iT{R -a-b))- fi2nT{R - a + b)) - f{2TxT{R + a - 6)) | , (4.6) 
where / is given by the power series (-B„ is the nth Bernoulli polynomial) 

°° 92n r> 

which satisfies the differential equation 

?/^/(l/) = cothy-^, /(O) = /'(O) = 0. (4.8) 
The differential equation may be solved numerically, yielding the results shown in Fig. [TTl 

V. NONCONTACT GEARS 

We consider first the interaction between corrugated planes, as shown in Fig. [121 Here 
we compute the lateral force Fiat between the corrugated plates, which are offset by a 



distan ce y n. T 
et al. 



33 



34 



le Dirichlet and electromagnetic cases were previously considered by Emig 



35|, to second order in corrugation amplitudes. We have carried out the 



calculations to fourth order 36|]. In weak coupling we can calculate to all orders, and verify 
that fourth order is very accurate, provided k^h <^ 1. We express the lateral force relative 
to the normal Casimir force between uncorrugated plates, F^^^: 

r = ^. ^ • (5.1) 

l-^Casl(^i^2/a^)/i;oa sm(A;oyo) 
The weak coupling limit is shown in Fig. [131 
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FIG. 11: Finite temperature interaction energy between two weakly-coupled semitransparent 
spheres, with radii a = b = R/4:, where R is the distance between the centers of the spheres. 
Shown are the exact result, the high T limit, and the truncated series expansion. Plotted is 
E = -XiX2abe/l6TrR. 




FIG. 12: Parallel, corrugated planes. 
A. Concentric corrugated cylinders 

We can also consider concentric corrugated cylinders, as shown in Fig. [T3] [31). For 
corrugations given by 5-function potentials with sinusoidal amplitudes: 

hi{9) = hismu{9 + 9o), (5.2a) 
h2{9) = h2smue, (5.2b) 

the torque to lowest order in the corrugations in strong coupling (Dirichlet limit) is 



r 

u sm z/i 



7T h h 1 
^O o^n ■^ ——B^J^^i'^)^ a = (as - ai)/(a2 + fli) i? = -(01 + 02). (5.3) 



2TrRL, "240a3 a a " ' ' " ' ' 2 
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FIG. 13: Weak coupling lateral forces for corrugated planes. Shown by the solid line is the exact 
lateral force, compared with the forces calculated in 2nd order and 4th order in the corrugation am- 
plitudes. Here we assume the two plates have equal corrugation amplitudes h. The oo superscripts 
signify that the result is exact both in h and a, relative to the wavelength of the corrugations. 




FIG. 14: Goncentric corrugated cylinders. 

Figure [15] shows the Dirichlet hmit of the torque on cyhndrical gears. A similar result can 
be found for weak coupling, which, again, has a closed form. 

We are currently completing analogous calculations for dielectric materials, illustrated by 
Fig. [T6l 2nd and 4th order results should appear soon, which will complement other recent 
work 
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FIG. 15: Torques between corrugated cylindrical gears calculated in 2nd order in the corrugation 
amplitudes, compared to the similar result for corrugated planes. 
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FIG. 16: Corrugated parallel dielectric slabs. 
VI. MULTILAYERED SURFACES 

An example of a simple multilayered potential is given in Fig. [171 To the left of this array 
of potentials, the reduced Green's function has the form, in terms of the reflection coefficient 
TZ for the array: 

g{x, x') = — ( e^"l"-"'l + 7^e"("+"')') , z,z' < 0. (6.1) 
2k \ / 

(We can actually find the Green's function everywhere, for any piecewise continuous poten- 
tial.) 

The array refiection coefficient may be readily expressed in terms of the reflection (R) 
and transmission (T) coefficients for a single potential: 

n = R + Te"""Re'""'{l - Re-'^'"'n)-\ (6.2) 

where a = \ — d is the distance between the potentials, each of thickness ci, and the result 
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FIG. 17: A semi-infinite array of periodic potentials. 



of summing multiple reflections is 



7^ 



1 

2R 



^2.a _l ^g2«a _ _ rpiy _ 4^2^-2 



(6.3) 



If the potentials consist of dielectric slabs, with dielectric constant e and thickness c?, the 
TE reflection and transmission coefficients for a single slab are (k' = a/^C^ + 

g2K'd _ ]^ 



R 



TE 



T 



TE 



A{K'/K)e 



K'd 



(6.4a) 



(6.4b) 



{i + K'/Kye^^'d-{i-K'/Ky 

The TM reflection and transmission coefficients are obtained by replacing, except in the 
exponents, k' n'/e. (Multilayer potentials have been discussed extensively in the past, 
see, for example, Refs. 28|, l39|, |40|, l4l|.) 



A. Casimir-Polder force 

Consider an atom, of polarizability <y{uj), a distance Z to the left of the array. The 
Casimir-Polder energy is 



E 



dc 



(27r): 



,a«)trg(Z,Z), 



where apart from an irrelevant constant the trace is 

tr g(Z, Z)^^ [-C'7^™ + (C' + 2^)71™] 



(6.5) 



(6.6) 
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FIG. 18: Casimir-Polder energy between an semi-infinite array of dielectric slabs with e = 2, 
compared to the energy (lower curve) if only one slab were present. Here we have assumed that 
the spacing between the slabs and the widths of the slabs are equal. 



FIG. 19: Casimir-Polder energy for large distances from the array, as a function of the ratio a/d, 
where a is the distance between the dielectric slabs in the array, and d is the thickness of each slab. 
Here e = 2. 

For example, in the static limit, where we disregard the frequency dependence of the polar- 
izability. 



This is compared with the single slab result in Fig. [181 It is interesting to consider the 
Z — > oo limit, which is shown in Fig. [191 When a/rf — > we recover the bulk limit. 

VII. ANNULAR PISTONS 

The multiple scattering approach allows us to calculate the torque between annular pis- 
tons, as illustrated in Fig. [201 We use multiple scattering in the angular coordinates, and 
an eigenvalue condition in the radial coordinates — the problem is equally well solvable with 
radial Green's functions, but this method permits interesting generalizations. 




10 



a /d 




(6.7) 
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FIG. 20: Two concentric Dirichlet cylinders, containing semitransparent radial planes separated 
by an angle a. The radial planes constitute an annular piston. 



b/a=: 

--- b/a=5 



FIG. 21: Energy of annular piston as a function of the angle between Dirichlet planes Ai = A2 = 00. 



Using the argument principle to determine the angular eigenvalues, we get the following 
expression for the energy between radial Dirichlet planes for an annular Casimir piston: 



00 



d 



Kdn / d7]—\n[Kin{Ka)Iin{Kb) - Iin{Ka)Kir^{Kb)] 

'i'^ i Jo J (^V 



7 



X In 1 



A1A2 cosh^ rj{TC — a)/ cosh^ rjir 



(7.1) 



{2ritanhriTC + Ai) (2?7 tanh 777? + A2 

where Itrj and Kij^ are modified Bessel functions of imaginary order, and the contour 7 
encircles the poles in rj along the positive real axis. As to be described elsewhere 42|, we 
have extracted numerical results for the energy, as shown in Fig. [2TJ 

We hope to apply these methods to study interactions between hyperbolae, such as a 
hyperbolic cylinder above a plane, and a hyperbola of revolution above a plane, see Figs 
and 
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FIG. 22: Hyperbolic cylinder above a plane. 



FIG. 23: Hyperboloid of revolution above a plane. 

VIII. CONCLUSIONS 

This survey describes some of the Oklahoma group's work applying multiple scattering 
techniques to interesting geometries. We have extensively explored weak coupling, because 
such cases can be carried out exactly, and therefore they serve as a laboratory for test- 
ing general features, such as edge effects and the range of validity of the proximity force 
approximation. We also have general results for the Green's functions for arbitrary piece- 
wise continuous potentials in separable coordinates. From these we can calculate not only 
Casimir-Polder forces, but Casimir energies and torques for many geometries, including 
annular pistons, and forces between hyperbolic surfaces. Finally, new results for electro- 
magnetic non-contact gears, both for conductors and dielectrics, are in progress. 
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